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\q . Abstract 

<n: 

We present a first numerical investigation of a non-commutative gauge theory denned via 
the spectral action for Moyal space with harmonic propagation. This action is approximated 
by finite matrices. Using Monte Carlo simulation we study various quantities such as the 



C3 



energy density, the specific heat density and some order parameters, varying the matrix size 
and the independent parameters of the model. We find a peak structure in the specific heat 
which might indicate possible phase transitions. However, there are mathematical arguments 
which show that the limit of infinite matrices can be quite different from the original spectral 
model. 

> 

in ■ 1 Introduction 
O 

cn . 

Quantum field theory on noncommutative spaces JlJ, [2|, [| is an active subject of research. The 
most-studied noncommutative spaces are the Moyal space [|J] and fuzzy spaces || . Fuzzy spaces 
are matrix approximations of manifolds and as such ideal for numerical investigations similar to 
non-perturbative quantum field theory on the lattice. In this paper we focus on the Moyal space, 
which is a continuous deformation of M. d for which the usual Fourier techniques of perturbative 
quantum field theory are available. It turned out that a renormalisable quantum field theory on 
M d is, in most cases, no longer renormalisable on d-dimensional Moyal space due to a phenomenon 
called ultraviolet/infrared mixing ||. In Q it was discovered that for the (^ 4 -model on 4- 
dimensional Moyal space the UV/IR-mixing generates an additional marginal coupling which 
corresponds to a harmonic oscillator potential for the free scalar field. The resulting action 

S[(p} = J d A x Q^*(-A + Sl 2 x 2 + fi 2 )-k(p + ^<p * <p * ip * <pj (x) (1) 

was then shown to be perturbatively renormalisable to all orders in A. In ([!]), x = -x, where 
G is the deformation matrix defining the Moyal product. See also |8|, ^, |l^, [ll], 12 1. Moreover, 



the frequency parameter can be restricted to Q 6 [0, 1] by Langmann-Szabo duality ]l3| 
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The renormalisability of the action (Q) raises the question whether a harmonic term can also 
render Yang-Mills theory renormalisable on Moyal space (recall that the usual Yang-Mills action 
on Moyal space has the same UV/IR- mixing problem JT3J). Yang-Mills theories in noncommut- 
ative geometry [15] are naturally obtained from the spectral action principle [16] relative to an 
appropriate Dirac operator. In |17]] it was shown that Moyal space (with usual Dirac operator) 
is a (non-compact) spectral triple; its corresponding spectral action was computed in Jl8| |. To 
obtain a gauge theory with harmonic oscillator potential via the spectral action principle, a 
differential square root of the harmonic oscillator Hamiltonian is necessary as Dirac operator. 



In absence of such a Dirac operator, in [19, |20[] an effective gauge model was constructed as the 
one-loop effective action of complex harmonic noncommutative quantum </? 4 -theory in a classical 
external gauge field. As a result, the noncommutative Yang-Mills Lagrangian is extended by 
two terms * and (X^ * X^) 2 , where X^ = ^x^ + is the 'covariant coordinate'. 

A first outline of a candidate spectral triple for harmonic oscillator Moyal space was given 
in [21]. Additionally, in [21] the linear and quadratic terms of the spectral action for a U(l)- 
Yang-Mills-Higgs model were computed and then extended by gauge invariance. Thereby the 
appearance of X^ * X^ was traced back to a deep entanglement of gauge and Higgs fields in a 
unified potential (aX^ * X^ + f3(p * (p — l) 2 , with a, /3 £ M + . 

It turned out that the candidate spectral triple proposed in |2l|] was the shadow of a new class 
of non-compact spectral triples with finite volume [22[|. The spectral geometry of Moyal space 



with harmonic propagation, which falls into this class, was fully worked out in 23]. There are 



in fact two (even, real) spectral triples (A*, H, T>,, T, J), with • £ {1,2}, for the c?-dimensional 
Moyal algebra A* and differential square roots "D. of the harmonic oscillator Hamiltonian. The 
spectral triples are of metric dimension d and KO-dimension 2d, have simple dimension spectrum 
consisting of the integers < d, and satisfy all regularity and compatibility requirements of spectral 



triples. Additionally, the spectral action was rigorously computed in [ 23 1 , i.e. with Holder type 
estimates for the remainder of the asymptotic expansion and with inclusion of the real structure 
J. 



A completely new feature of the spectral action [23, (and also of the effective action 
[19, |2(J) is that the expansion of X^ * X^ and its square produces a term which is linear in 



the gauge field A. This means that the vacuum, i.e. the solution of the classical field equations, 
is no longer taken at A^ = (or more generally at a flat connection Fa U = 0) but at some 
non-constant value for the gauge field. A first discussion of the vacuum structure of this type of 
gauge models was given in [p4| . It turned out that generically there are infinitely many vacuum 
solutions. Some of them were exposed, but it was not possible to give reasonable argument for 
the right solution. In particular, it became completely impossible to study the gauge model as 
a perturbative quantum field theory. 

This is the point where the numerical treatment comes into play. The standard method 
of numerical quantum field theory is to approximate the space by discrete points, for example 



using a lattice approximation and then calculate the observables over that set of points [25]. For 
Moyal space a position space approximation is not suitable due to the oscillator factor of the 
Moyal product. Instead, we shall use the matrix Moyal base (which was already used in the first 
renormalisation proof j?]] of (/? 4 -model), restrict it to finite matrices and perform a Monte Carlo 
simulation of the resulting action. In this way we will study some statistical quantities such as 
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energy density and specific heat, varying the parameters Q, ^-r-,a of the model and gathering 

xo 



X-i 

xo 

some information on the various contributions of the fields to the action. The simulations are 
quite cumbersome due the complexity of the action and the number of independent matrices 
to handle. Nevertheless we are able to get an acceptable balance between the computation 
precision and the computation time. For the simulations we apply a standard Metropolis-Monte 



Carlo algorithm [26] with various estimators for the error and for the autocorrelation time of the 
samples. The range of parameters is chosen to avoid problems with the thermalisation process, 
thus permitting a relative small number of Monte Carlo steps to compute independent results 
from the initial conditions. 

We are eventually interested in the continuum limit which corresponds to matrices of infinite 
size. We thus compute our observables such as the energy density for various matrix sizes and 
then look for a stabilization of these observables as the matrix size increases. The specific 
heat, which is a measure of the dispersion of the energy, will be used to identify possible phase 
transitions in form of peaks of the specific heat at increasing matrix size. 



2 Four-dimensional harmonic Yang-Mills model 

The harmonic Yang-Mills model is defined as the spectral action resulting from the spectral 
triples (A*,7-L,U,,r, J), with • G {1,2}, analysed in |23|] . The Moyal algebra .4* is the space of 
Schwartz class functions on M 4 equipped with the product 



f*g(x)=[ 7t4t f{x+\Q ■ k)g(x+y)e 



Where {k,y) is the scalar product defined as k^y^. The unbounded selfadjoint operators T>, 
on the Hilbert space % are differential square roots of the harmonic oscillator Hamiltonian 
H = —d^d 11 + frcc^x^ of frequency Cl, i.e. T> 2 = H — (—l)°£lT>, for a certain spin matrix S. If 
L*(f) denotes left Moyal multiplication with a function / G A±, then one has 

[Z>i, £*(/)] = LSd^)®^ , [V 2 , £*(/)] = L,(^/)®r^+ 4 , (3) 

where the matrices Ti, . . . , Tg satisfy the anticommutation relations 

{i* r^} = {r^+ 4 , r u+4 } = 2{ g ~ 1 Y u , {r^, r^ 4 } = o , (4) 

relative to an induced metric g = (id — \VL 2 Q 2 )~ l . The grading is V = ri-'-Tg, and the 
real structure satisfy JP,J _1 = V, and JL ir (f)J~ 1 = R±(f), where -R* denotes right Moyal 
multiplication. 

In order to implement the Higgs mechanism a la Connes-Lott 27] one considers the product 



of the spectral triple (A*, H, T>i,T, J) with the finite Higgs spectral triple (C © C, C 2 , Mai, Jf), 
where o\ is a Pauli matrix, Jf any matricial real structure and M > 0. Then, a self-adjoint 
fluctuation A = ^2ai[V,bi] of the total Dirac operator V = {V\ © 1 + V © Ma\) to give 
T> A = V + A + JAJ" 1 , for J = J © J/, is of the form 
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for the components A^,B^ E A.* of two real one- forms and a complex scalar (ft S A+. Using T>2 
instead of T>\ amounts to replace by T^ +4 . 

The spectral action principle [16] asserts that the bosonic action of a field theory with 
fermionic Dirac operator T>a has the form 

S(V A ) = Tr( X (V A )) , (6) 

where x is a smooth approximation of the characteristic function on [0, A 2 ], for some scale 
parameter A. For the fluctuation (||), the part of the spectral action which is relevant and 
marginal for A — > oo has been explicitly computed in |23], for general effective metric g. This 
computation involved Laplace transformation, Duhamel expansion with Holder-type estimates 
for the remainder and explicit use of the Mehler kernel for the harmonic oscillator Hamiltonian. 
For a special choice of the noncommutativity matrix 2 = — # 2 id, the result of [23| takes in 
terms of := ^p, the moments A 2n x~n := J °° ds s n ~ 1 x(s) of the "characteristic function" and 
Xo := x(0) the form 

M 2 6> 4 A 6 / M 4 fl 4 A 4 80 2 A 4 \ _ / M 6 6 )4 A 2 2M 2 6> 2 A 2 \ 

-52 M 8 6> 4 M 4 fl 2 \ 
,45 + 192iF + ~3O 2_ ; X0 

vr 2 (l + 2 ) 2 / ^(l + O 2 )^*^ 

/ AO 2 a/ \ 2 / AV) 2 -v 

+ + M(0+0) + -^f-X^*^ + M 2 — *^A 2 ) - (-^f_Xo^ + M 2 — *^A< 
V 1+U Z A xo ' vi+sz 2 ^ xo 

4^ 2 ^ , „,2 X-i a2 \ 2 / 40 2 - ~ H , , „ x-l 



+ Xo 



+ ( + M(cP+4>) + ^—X B *X* + M 2 — ^A 2 - __X *x£ + M 2 - ^A< 

v i+n 2 p xo y vi+o 2 p xo 

+ 0(A~ 1 ). (7) 

Here, = d^cj) — iA^ * <j) + i(f> * -B^ — iM(A fl — B^) is the covariant derivative of the scalar 
field, := d^A v — d v A^ — i{A^ * A u — A v * ^4^) the field strength of ^4 and similarly 
the field strength of B. Moreover, Xau '■= Xqu + A^ and Xbu '■= Xqu + B^ are the covariant 
derivatives of A and 5, respectively, where X ^ ■= ^- = (0 _1 ) ^x" '. The remarkable outcome 
of the spectral action (|7|) is that the Higgs field <f> and the gauge fields A, B appear together 
in a unified potential. In this way, also the gauge field shows a non-trivial vacuum structure. 
Besides, the action is invariant under U(A-k) x U(A*) transformations: 

{(f) + M) ua* {4> + M) * Tig, X^ ^ua + X^ + ua, X^ B ^ u B -kX^*W§ . (8) 



3 Discretisation by Moyal base 

The 2-dimensional Moyal algebra with deformation parameter 9 > has a natural basis of 
eigenfunctions f mn of the harmonic oscillator, where m,n £ N. These are given in radial 
coordinates by 

f mn (pcoB<p,psm<p) = 2(-ir^e^"-™) U^pj e^L^™ (y 2 ) (9) 
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and satisfy 

$nkfml(x) (10) 

2ir95 mn , (11) 

for details. The expansion of Schwartz functions on R 4 in the Moyal base, 
A* 3 a = a(x , ■ ■ ■ ,x 3 ) = } a m i n i fmmi(x ,xi)f m2n2 (x2,x 3 ) , (12) 

< * m 2 n 2 
mi,m2,ni,n26N 

then provides an isomorphism of Frechet spaces between A* and the space of rapidly decreasing 
double sequences (a mn ) m ngN 2 equipped with the family of seminorms 

OO 1 

Pk((a mn ) m ,neN 2 ) ■= (( 2 H + 1 ) 2 ' C ( 2 l n l + l) 2fc |a»nn| 2 ) 2 , \m\ := m 1 + m 2 ■ (13) 

m,n£N 2 

According to dlO|), Moyal product and integral reduce in the (/ mn )-basis to product and trace of 
infinite N 2 -labeled matrices, with convergent index sums due to (|l3|). By duality, the covariant 
derivatives and can also be expanded in the (/ mn )-basis, but the expansion coefficients 
X^ mini , X^ mini diverge for m,i,ni — > oo. 

To any a S A* we can associate a sequence (cl n )n£N of cut-off matrices 

if max(mi, r/12, n>i, n 2 ) < N , 
else . 

Then, (a N ) is a Cauchy sequence in any of the semi- norms pk and converges to a in the Frechet 
topology of A*. 

In quantum field theory we are confronted with the converse problem. To deal with diver- 
gences, a regularisation has to be introduced which restricts the system to a finite number of 
degrees of freedom. After re-normalisation from bare to physical quantities one has to show that 
the limit to an infinite number of degrees of freedom is well-defined. In our case, the natural 
regularisation is to restrict the matrix indices to m; < N, which corresponds to a cut-off in the 
energy. Even if we could solve the renormalisation problem, the removal of the cut-off, i.e. the 
limit N — > 00 to infinite matrices, will fail: A sequence of (N x iV)-matrix algebras does not 
converge in the Frechet topology. 

Fortunately, in quantum field theory we are interested in the convergence of correlation 
functions, and not of matrix algebras. The path integral in usual quantum field theories is over 
random walks and not over smooth field configurations. It seems not impossible (although we 
cannot prove it) that quantum correlation functions are less sensitive to the topology of the 
underlying classical field theory. In our case, we make the (not verifiable) hypothesis that the 
cut-off correlation functions to be computed carry some information about the original smooth 
model. 

Using the identities = i{4>+M)-kX B ^-iX A ^((j)+M) and F£ v = -i[X$ , X^]+i[X^, X%] 
(and similarly for F^), and ignoring all contributions of X® which for finite matrices yields some 



(fmn*fkl)(%) = 

J d %fmn(%) — 



N 

m 2 n 2 



m 2 n 2 
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finite number, we can recast the restriction of the action (^) to finite matrices in the following 
form: 



s{4>,x A ,x B ) 



1 



i_o 2 ) 2 (i + n 2 ) 



+ 



(1 + ft 2 ) 2 



Tr<( 



6(1 + ft 2 ) 2 



X AfJi , Xau 



+ 4>-k^ + M{4>+4>) + 



4 ^\j^*X^ + M 2 -A 2 *^ 

Xo 



1 + S1 2 A 



+ 4>*<j) + M \<ft+<f>) + 



aq 2 



Xn -k Xbu + M — A 



1 + fi 2 s 



+ 2(1 + n 2 ) ( (0 + M) * x Bii - x All * {cb + M) 



2 X-1 

Xo 



((0 + M) * - X£* (0 + M) 



(14) 



The restriction to finite matrices shows crucial differences to the smooth model. Only these 
differences make the numerical simulations possible with the drawback of the serious possiblility 
that our results can deviate from the original smooth model. 



1. The action (14) has an obvious family of minima given appropriate multiples of the identity 
matrices. We thus define 



+ M = ip + A A /— cos al 

Xo 



2Q 2 



X A » = Y A » + \ A^W i M siD n 



X B u = Y B u + -A 



1* fx 



2Q 2 



2 V Xo V( 1 + ^ 2 ) " 



I,, sin a . 



(15) 
(16) 
(17) 



Note that the corresponding minimum configurations for A^,B^ explicitly violate, in the 
limit JV — > 00, the Frechet condition. 

Substituting the previous fields into fll4| ) we get a positive action with minimum in zero: 



S(ib,Y A ,Y B ) 



1 



Tr D[ [Y Aix ,Y Av ] [YX,YX] + [Y B „Y Bu ] [Yg,Yg] 



(i + n 2 ) 2 

+ (ibtp + [i cos a(ib + + CY^Yavl + Ml"^ sin a) 
+ (^V + A* cos «0 + VO + CY^Y Bll + /xI m Yb m sin a) 

+2(1 + n 2 ) ((Yb m - r^j/i cos a + ^y S/i - y a ^) 

((Y£ - Yg)fi cos a + $YX- Ygj) \ , 



(18) 



with 



C 



1 + n 2 
~lv 2 ~ 



D 



'1-n 2 ) 2 (1 + 2 ) 4 

2 6(1 + ft 2 ) 2 



A : 



2 X-1 
Xo 



(19) 
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2. For finite matrices, the N 2 -indexed double sequences can be written as tensor products of 
ordinary matrices, 



K 



x™ini = vr mini ®X % . (20) 

mono ' ' 44 



m 2 n 2 

1=1 



Since the matrix product and trace also factor into these independent components, the 
action factors into S = Y,?=i S(ip u , Y% , Y^)S{ip 2i , Yf, Yf ). Then, regarding all ip u , Y£\ 
Y B , ip 21 , Y A , Y b as random variables over which to integrate in the partition function, the 
partition function factors, too: 



/-TV /ll vll vll /21 v-21 v2l\ T<( ilK V 1K V 1K ,2K V 2K V 2K\ -S 

/ V{V ,Y A ,Y B ,ip , Y A , Y B ) • • • P(V> , Ya , Y B , ip ,Y A ,Y B ) e 
= ^jv{^\ Y\\ Yl\i?\ Yf\ Yf ) e - s ^ U ' Y A' Y B )-S(^,Y»Yg) 



(21) 

We may therefore restrict ourselves to K = 1. Now the discretized action is invariant 
under the same transformations (||), but now with (u A ,u B ) £ U(M.jy) x f/(Mjv). For the 
limit JV->oo we would have K — > oo and therefore problems with convergence. 

3. Instead of integrating in the partition function over all gauge-equivalence classes of ip, Y A , Y B 
as required, we follow the usual matrix model philosophy^ and integrate over all matrices 
ip,Y A , Y B . A reduction of this 10-matrix model to gauge-equivalence classes seems rather 
hopeless. 



It is convenient to pass, for each factor in the tensor product (20), to complex matrices 



Z 


= Y A + l Y 1 A , 


z 


- Y A 


-*Y A 




= Y Q B + iY*, 


Zi 


- Y B 


-iYf 


z 2 


= Y 2 A + iY 2 A , 


Zi 


- Y A 

— r 2 


— iY A 


z 3 


= Y 2 B + iY z B , 


z 3 


- Y B 

— r 2 


~ *Y* 



(22) 

The convention that the bar denotes the hermitian conjugate will also be used for the complex 
matrix ip. In the end using the substitutions ( ^2|) and after some simple manipulations, the 
discretised action is: 

5*4 = Tr [Cf + Cvo + £vi + £d £-Dq + £d 1 £d 1 + £d 2 £-d 2 + £d 3 Ad 3 ) , (23) 



with 



C F = ^r{ [Z , Z ] 2 + [Z x , Z{f + U [Z + Z , Z 2 -Z 2 ] 2 - [Z + Z , Z 2 + Z 



2 y L u ' U J la' x J 4 

+ [Zq — Zo, Z 2 + Z 2 ] — [Zq — Zo, Z 2 — Z 2 ] — \Z\ + Z\,Z 3 + Zs\ 
+ [Z x + Z x , z 3 -z 3 } 2 + [Z x - Zi, z 3 + z 3 } 2 - [Zi - Z u z 3 - z 3 ] 2 



2 

21 



1 In scalar 1-matrix models the gauge-equivalence classes are the configurations of eigenvalues. Nevertheless it is 
custom to integrate over all matrices. The integration over the gauge group produces a measure for the eigenvalues 
which is given by the square of the Vandermonde determinant. In 2-matrix models the gauge equivalence classes 
are the eivenvalues together with a unitary matrix which describes the relative orientation of the eigenbases. 
The full matrix integration can be reduced to an integration over the eingenvalues thanks to the formula of 
Itzykson-Zuber and Harish- Chandra. 
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£v = {ipip + V cos a(ip + + - ({Z ,Z } + {Z 2 ,Z 2 }) 
+ ^^((-1 + i)(Z + Z 2 ) + (1 + i)(Z + Z 2 ))f 
C Vl = (ipip + /j,cosa(tp + i>) + 2 {{Zi,Zi} + {Z3,Z 3 }) 

+ ^^((-1 + 0(^1 + 3») + (1 + + ^3))) 2 

£ Do = ^2(l + 2 ) (// cos a(Zi +Z 1 -Z - Z ) + +Zi)-(Z + Z )if>) 
C Dl = y/2(l + ffl) (/i cos a(Zi - - Z + Z„) + - Z x ) - (Z - 
£ fl2 = V / 2(TTfF) (/i cos a(Z 3 + Z 3 - Z 2 - Z 2 ) + ^j(Z 3 + Z 3 ) - (Z 2 + Z 2 )^) 
£z3 3 = y/2(l + Q 2 ) (n cos a(Z 3 - Z x - Z 2 + Z 2 ) + V(^3 - 30 - (Z 2 - Z 2 )^) 



In this case, (23) becomes an action for 5 complex matrices. For the partition function we 
need the independent product of two copies of ([23]), i.e. we are dealing with a complex 10-matrix 
model. This is already cumbersome and shows that there is little hope to treat the original model 
(0) where the simplifying consequences of finite matrices are not available. 

The next step is to define the estimator for the average values of interest and to specify some 
numerical parameters in order to analyse the numerical results. 



4 Definition of the observables 



Following Monte Carlo methods, we will produce a sequence of configurations {(ip, Zi)j}j = i i2t ... ^ MC 
and evaluate the average of the observables over that set of configurations. These sequences of 
configurations, called Monte Carlo chain, are representatives of the configuration space at given 
parameters. In this framework the expectation value is approximated as 

<°> « r— E °i . ( 24 ) 

where Oj is the value of the observable O evaluated in the j-sampled configuration, (if),Zi)j, 
Oj = 0[(ip, Zj)j\. The internal energy is defined as 

E(n,f,,a) = (S) , (25) 

and the specific heat takes the form 

C(n,»,a) = (S 2 )-(S) 2 . (26) 

These quantities correspond to the usual definitions for energy 

1 87. 

m»,<*) = ~^ (27) 

and specific heat 

8 W 

C(n,»,a) = — , (28) 
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where Z is the partition function. It is very useful to compute separately the average values of 
the four contributions: 



F(n,fi,a) = (TrC F ) , (29) 

V (n,ti,a) = {TrC Vo ), (30) 

Vi(n,fj,,a) = (TrC Vl } , (31) 

D(Q,n,a) = (Tr [C Do C Do + • • • + £ Ds £d 3 ) ) • (32) 

4.1 Order parameters 

The previous quantities are not enough if we want to measure the various contributions of 
different modes of the fields to the configuration (ip,Zi). Therefore, we need some control 
parameters usually called order parameters. As a first idea we can think about a quantity 
related to the norms of the fields, for example the sums ^2 nm (V'nmpj ^2inm \Zinm\ 2 • These 
quantities are called the full-power-of-the-field [^, ^9| ; they can be computed as the trace of 
the square: 

= Tr(|Vf) (33) 
Z^ = Tr(|^| 2 ) (34) 

In contrast, (tp a ) alone is not a good order parameter because it does not distinguish contribu- 
tions from the different modes. But we can use it as a reference to define the quantities 

N 

2 _ I |2 
fO ~ / j \ a nn\ > 

n=0 

N 

%i0 = ^ ] \ z inn\ , (35) 
n=0 

where a mn and Zi mn are the expansion coefficients of ip and respectively, in the matrix base 
@ . Referring to @ it is easy to see that these parameters (|35|) are connected with the purely 
spherical contribution. These quantities will be used to analyse the spherical contribution to 
the full-power-of-the-field. We can generalise the previous quantity and define parameters (pi in 
such a way that they form a decomposition of the full-power-of-the-fields: 

Z>0 Z>0 
Following this prescription, the other quantities for / > can be defined as: 

l l 

<Pl = ^ J | Oram (1 _ <W)| 2 , %il= X/ \ z lnm{^ - 5 n m)\ 2 ■ (37) 

n,m=0 n,m=0 

If the contribution is dominated by the spherically symmetric parameter we expect to have 

In the next simulations we will evaluate, apart from I = 0, the quantity with / = 1 as 
representatives of those contributions where the rotational symmetry is broken. According to 



(|37D we have 

fi = \ a io\ 2 + |aoi| 2 , Zfi = kiiol 2 + \zm\ 2 ■ (38) 
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Using higher I in ( |37| ) we could analyse the contributions of the remaining modes, but it turns 
out that the measurements of the first two modes are enough to characterise the behaviour of 
the system. 

5 Numerical results 

Now we discuss the results of the Monte Carlo simulation of the approximated spectral model. As 
a first approach we use some restrictions on the parameters. Starting point is the approximation 
( p3|) of the spectral action. Since ( p3| ) is symmetric under the transformation \i i— >■ — \i we can 
assume \i > and /j, 2 > 0. In this first treatment we explore the range [i G [0,3.1], which is 
enough to show a particular behaviour of the system for fixed 0. The parameter appears only 
with its square and is defined as a real parameter, therefore also for Q we require > 0. For the 
scalar model @ it was possible to restrict to G [0, 1], because Langmann-Szabo duality maps 
Q to jj. In the gauge model under consideration, Langmann-Szabo duality is not realised. Due 
to the prefactor in front of the integral (p3|), the action vanishes for $7 — > oo. Studying the plots 
for the energy and the specific heat, we have chosen the range f2 G [0, 2tt] in which the action 
is significantly different from zero. The last parameter to consider is a, which is connected to 
the choice of the vacuum state, with range a G [0,2-7r]. The study of the system varying a is 
quite important from a theoretical point of view because it is related to the vacuum invariance. 
In the action there appear some contributions proportional to (sin a)/Q, which seem to diverge 
for rZ = 0. Numerically we have verified that this is an eliminable divergence and the curves 
of the observables can be extended to f2 = by continuity. Studying the dependence on a we 
can conclude that in the limit N — > oo the observables are independent from a, therefore for 
our purposes a will be fixed equal to zero avoiding the annoying terms. In general, for each 
observable we compute the plots for matrix size approximations ./V =5, 10, 15, 20. 

5.1 Varying a 

We start looking at the variation of the energy density and of the full-power-of-the-fields density 
for fixed \i and Q, varying a G [0, 2tt}. As representatives we present the plots for [i = 1, 
Q G {1, 0.5}, but we obtain the same behaviour for any other choice of the parameters allowed 
in the considered range. All three plots show an oscillating behaviour of the values, and this 




.1 4 5 6 113 4 5ft" I 2 .1 4 5 



Figure 1: Total energy density and full-power-of-the-fields density for {(fa), (ZQ a ) (from the left to the right) 
fixing = 1, Q = 1, varying a and N. N = 5 (circle), N = 10 (square), N — 15 (triangle), TV = 20 (cross). 

oscillation is present in all other quantities measured. The amplitude of this oscillation becomes 
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smaller and smaller increasing the size of the matrix and this is true for all the quantities. The 
same trend is described in fig.|| which shows different positions of the maxima, but again smaller 
amplitudes for increasing N. These results allow us to consider a = for all next plots, since we 




f, " 1 2 .1 4 ? 6 " 12 3 4 5 6 



Figure 2: Total energy density and full-power-of-the-fields density for (ifia), (Zo a ) (from the left to the right) 
fixing = 1, Q = 0.5, varying a and N. 

are interested in the behaviour of the system for N — > oo. This occurrence simplify all the next 
simulations thanks to the vanishing of terms ~ (sin a) /Vt appearing in the discretised action. 

5.2 Varying Vt 

As already mentioned we chose [0, 3] as range for Q,. In fact, if we look at the plots in figH of 
the total energy density (S)/N 2 for \i £ {0, 1}, we notice that the action tends to zero for Q 
outside the selected interval. This behaviour of the action is the same for all possible choices of 
parameters and for the specific heat, too. 



El N- EIN 1 




n - 

0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 3: Total energy density and the various contributions for /i = (left), fi — (right), a = varying Q 
and N. With N = 5 (circle), N = 10 (square), N = 15 (triangle), N = 20 (cross). 

In the rest of this section we ignore for the computations of (E), (D), {V}, {F) the global 
prefactor (1 + f2 2 ) -1 . In this way we focus our attention to the integral as the source of possible 
phase transitions. Now we will analyse three cases in which \x is fixed to 0, 1,3. In all cases a 
is zero and we vary Q £ [0, 3]. The plots in fig.§ show the total energy density and the various 
contributions: the potential V/N 2 , the Yang-Mills part F/N 2 and the covariant derivative part 
D/N 2 , for fj, = 1. There is no evident discontinuity or peak, and increasing the size of the 
matrices the curves remain smooth. 

Comparing the energy density and the various contributions in fig.[| we notice that the 
contributions between F and V balance each other and the total energy follows the slope of D, 
and this behaviour continues increasing the size of the matrices. 
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Figure 4: Total energy density and the various contributions for \i = 1, a = varying Q and N. From the left 
to the right E, V, D, F with N = 5 (circle), N = 10 (square), AT = 15 (triangle), JV = 20 (cross). 



Ent:/,V- 




0.5 1.0 



2.0 2.5 3.0 



0.5 1.0 1.5 2.0 



Figure 5: Comparison of the total energy density and the various contributions for = 1, a = 0. E (circle), F 
(triangle), D (cross), V (square). With N = 5 (left) and N = 20 (right). 



The specific heat density in fig.^ shows a small peak in Q = 0. This peak does not increases 
as TV increases, therefore is not related to a phase transition. 




2.0 2.5 3.0 



Figure 6: Specific heat for fi = 1. 



In order to gain some information on the composition of the fields we look at the order 
parameters defined in the previous section. Starting from the scalar field ip, fig-^ shows the 
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plots for ((fa), (ifo) and (ipf) for N = 5. The three values {^>q) and (iff) seem essentially 

constant, where the spherical contribution (<^>q) to the full-power-of-the-field is dominant. The 




Fi gure 7: On the left comparison of (tp\) (circle), (9^0) (square) and (y?i) (triangle) density. On the right 
comparison of (Zo a ) (circle), (Zoo) (square) and (Zqi) (triangle) density. 

behaviour of the Zq fields (which describe the covariuant coordinates) is different. Here the 
spherical contribution becomes dominant only for fi approaching 0, starting from a zone in 
which the contribution of (Zq ) and (Z^) are comparable. For brevity we only show the plots 
for {ZQ a ), (Zq ) and (Z^), but taking into account the statistical errors, the other Zj-related 
plots are compatible to the Zo-case. The dependence of the previous quantities on N is shown 
in the following plots fig§ All previous parameters decrease with N, but the dominance of (po 



<3>IH- <¥5>/N 2 «fi>IN 2 




Fi gure 8: Starting from the up left corner and from the left to the right the densities for {<fia)i (Po)i Z$ a , 
{Zqo} an d (Zoi) for fi = 1 varying SI and N. 

on the total-power-of-the-field is independent by N. The peak related to Zq decreases with N, 
but if we look at the single plot for the spherical contribution at N = 20, the peak persists as fi 
approaches fi = 0. 

Now we will analyse the model for /i = 0. Fig]9| shows the plots for total energy density and 
the contributions V, D, F. The slope of the total energy density seems to be constant. The 
^-contribution and the .F-contribution do not balance each other like in the previous case, but 
all three contributions balance themselves to produce a constant sum. 
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0.5 1.0 



Fi gure 9: Total energy density, various contributions and the comparison among them for — varying SX and 
N. From the left to the right E, V, D, F and comparison. 



The specific heat density fig. 10 shows again the small peak in Q = without iV-dependence. 




1.0 1.5 2.0 2.5 3.0 



Fi gure 10: Specific heat density for \i = varying £1 and N . 



For the other quantities (</?o)> (Vi) an d (-^oa); (^oo)> (^oi) we have according to fjg.pTl] 

the same behaviour as in the case [i = 1. 

A completely different response of the system is obtained in the plots for /i = 3, as we can 
see from fig.|l2| The slope of total energy density is very similar to the .F-component instead of 
D. However, there appears a sharp minimum around Q = 0.1 and two maxima at O ~ 0.6 and 
~ 1.8 for large N. This dramatic change in the plots might be interpreted as consequence of 
a phase transition in the parameter /u. Actually, in the next section we will find a peak in the 
specific heat density for some fixed £1 and varying [0,3]. 



The specific heat density fig.|13| displays a strong change, too. In fact, instead of the peak 
at $7 = 0, the peak appears close to the origin around Q = 0.15. This peak, in contrast to the 
previous ones, grows as N increases and therefore could indicate a phase transition. 

The fig. 14 describes the behaviour of the order parameters densities (<^), (^>q), and 
(Zq q ), {Zqq), They show a similar behaviour as the corresponding plots for fi = 1 and 

fi = 0. For the ip field the spherical contribution remains dominant. However, in the {(p\) plot 
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Figure 11: Starting from the up left corner and from the left to the right the densities for {ipa}, (<Po), {<Pi}, 
(Zoa), (Zoo) an( i {Zqi} for = varying fl and N. 




n 



0.5 ID 15 2.0 2.5 3.0 0J 1.0 1.5 2.0 2.5 3.0 



Figure 12: Total energy density, various contributions and the comparison among them for [i — 3 varying SI 
and N. From the left to the right E, V, D, F and comparison. 

there appears a deviation from the constant slope. This deviation is evident for N = 5 but still 
present for higher N. The order parameters for Zq display a peak close to the origin without 
oscillations even for N = 5. This maximum for higher N does not move closer to the origin, in 
other words, this shift is not caused by finite volume effects. Even for there appears a peak 
at f2 = which becomes shifted and smoother for higher N. 



5.3 Varying fi 

In this section we analyse the response of the system varying /j, £ [0, 3] while f2 is fixed at 0, 1 or 
3, and a is always zero. We start displaying the plots fig. 15 of the total energy density and of 
various contributions for Q = 0. There is no evident discontinuity but there appears a peak in 
the total energy density around [i ~ 2.5 for ./V = 20. Comparing all the contributions it is easy 
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Fi gure 13: Specific heat density for \i = 3 varying £1 and N . 



to notice that the slope of the total energy is dictated by the curve V of the potential part. 

As mentioned before, the specific heat density fig. 16 features a peak around fi ~ 2.5 for 
N = 20. Again, since the peak increases with N, we could relate this to a phase transition. The 
plots for the quantities (y>„) and (</?q) show a strong dependence on /u, in particular the slope of 
(c/Jq) seems mostly linear. The plot for ((pf) also increases with /u, but not linearly. From the 
first three plots of fig.0 we deduce that close to the origin the non-spherical contribution (ipf) 
is bigger than the spherical one (</2q). Increasing fj,, this situation capsizes and (tpo) becomes 
dominant over (fi). 

The behaviour of the Zq fields as shown in the last three plots of fig.|l7] is quite different. 
The spherical contribution is always dominant for the whole interval \i 6 [0,3]. The curves for 
(ZQ a ), (Zq ) are compatible to the constant slope. For (Z^) we have the same dependence on fj,, 
in particular there is a smooth descending step which becomes smoother for bigger N, However, 
we admit that due to some cancellation effects, the statistical errors are quite big so that this 
interpretation is not fully conclusive. Anyway, this result demonstrates the dependence of the 
order parameter for Z{, and in general of the system, on the two choices 17 = or ^ 0. 

Now we will analyse the model for Q = 1. As fig. [IB shows, the plots have a different slope 




05 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 25 3.0 



Figure 14: Starting from the up left corner and from the left to the right the densities for (<Pa), (<Po), (<Pi), 
{Zoa}, { z oo) an d { z oi) f° r /" = 3 varying Q. and N. 
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Figure 15: The total energy density and the various contributions for $1 = varying /i and TV. From the left to 
the right E, V, D, F an comparison with TV = 5 (circle), TV — 10 (square), TV = 15 (triangle), TV = 20 (cross). 
For the comparison: E (circle), V (square), D (triangle), F (cross). 

C/N 2 




0.5 1.0 1.5 2.0 2.5 3.0 



Fi gure 16: Specific heat density for $7 = varying and TV. 



compared to the previous case. The maximum of total energy density follows the one of the 
V^-component. If we focus only on the total energy plot and compare it with the one for O = 0, 
we notice a shift of the maximum for each TV. In particular, in fig.pl^ some maxima are moved 
outside the considered interval. We can find this shift very clearly looking at specific heat density 
plotted in fig. 19. Here again the peak both increases with TV and is shifted to \i ~ 3.3. 

Fig]2^ shows for (ifo) the same behaviour as in the case Q = 0. The plot for (<pf) 

displays an almost constant curve. However, close to the origin, the spherical contribution and 
the first non-spherical one are comparable. The introduction of f2 ^ creates, in the Zo-order 
parameters shown in fig.^O], a dependence similar to the plots for ip. The full-power-of-the-field 
density and the spherical contribution are no longer constant, they grow as [i increases. Even 
in this case the spherical contribution is always dominant excluding the region around /x = 0. 

The last set of plots treats the case = 3. The following diagrams for the energy and its 
contributions show the absence of the previous peak. They show a sort of dilatation of the 
former plots of fig,|T§|, 



The specific heat density fig. 22 does not show the peak in zero anymore, and the curves do 
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.0 1.5 2.0 25 3.0 



0.5 1.0 1.5 2.0 25 iD 



Figure 17: Starting from the up left corner and from the left to the right the densities for ((fa), (ipo), (fi), 
{Zl a ), (Zoo) an d ( z ai) Ior ^ = varying fi and N. 




.0 15 2.0 25 3.0 



Fi glire 18; Total energy density and contributions for Q = 1 varying fi and N . From the left to the right V , 
D, F. 




1.0 1.5 2.0 2.5 3.0 



Fi gure 19: Specific heat density for £1 = 1 varying fi and iV. 
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Fi gure 21; Total energy density and the various contributions for £1 = 3 varying ^ and N . From the left to the 
right E, V, D,F . 

not show any particular point as N increases. Actually, the peak can be found for higher fi. 

At last in fig.^3| we find a behaviour of the density of the order parameters for the Zq and 
ip fields similar to the former plots for = 1, and they are compatible with a dilatation of the 
previous plots. 



Conclusions and prospectives 

We have studied a noncommutative gauge theory which arises by restriction of the spectral 
action for harmonic Moyal space to finite matrices. For this quantum field theoretical model 
we have performed Monte Carlo simulations to obtain, as function of the parameters (//, Jl), 
non-perturbative information for the energy density, for various contributions to the energy 
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Figure 23: Starting from the up left corner and from the left to the right the densities for (ipa), (vo)) (v 9 ?); 
{Zoa}, (Zoo) and (Zoi) f° r ^ = 3 varying \i and N. 

density and for the specific heat density, as well as for a set of order parameters related to 
sphericity. Despite the complexity of the approximated spectral action considered here, we were 
able to obtain some reliable numerical results, showing that a numerical treatment of this kind 
of noncommutative gauge models seems feasible. However, as the restriction to finite matrices 
shows severe differences to the original smooth action, the relevance of our results to the smooth 
case is not clear. 

The specific heat density shows various peaks which could indicate phase transitions. In 
particular, studying the behaviour for some fixed mass parameter fx we found a relevant peak 
close to f2 = for fx = 3, and we noticed a big change in the energy density and in its contri- 
butions between the cases \x £ {0, 1} and fx = 3. Other peaks in the specific heat density can 
be found varying fx and fixing 0,. The plots show that increasing £1, the peak in the specific 
heat which starts at fx ~ 2.4 for f2 = is moved towards higher fx. The order parameters we 
introduced show a strong dependence on the cases £1 = versus f2 7^ 0. Referring to the fixed-/i 
plots we found a peak in the spherical contribution for the gauge fields Zj. Its behaviour can be 
interpreted as a sort of symmetry breaking introduced by 0, 7^ 0. Additionally, varying fx and 
fixing $7, the other parameters display a slope increasing with fx for all fields and all situations 
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but one: the plots of the order parameters (Zo a ), (Zqq) for O = show a constant behaviour. 

The natural next steps in the numerical study of this model could be the computation of 
the transition curves in order to separate the phase regions and to classify them using possibly 
additional order parameters. Our treatment, forced by limited resource, was conducted varying 
f2 in the range [0,3.1], since the Langmann-Szabo duality does not hold anymore in our case. 
Actually, the computed plots do not show any periodicity in f2 G [0, 1] so that we can infer that 
in contrast to the scalar case the range [0, 1] is not enough to describe the system. 

It will be very interesting to relax the condition /i 2 > 0. Implementing /i 2 < amounts to 
conduct the calculation no longer around the minimum of the action, in particular avoiding the 
explicit use of the finite vacuum and considering directly the four indexed components we can 
hope in some improvement in the behavior of the large N limit making the simulations much 
more harder (but not impossible), however the limit can continue to show big differences. Such 
continuous limit issues are not new for the simulations of a scalar field theory on matrix model, 
in facts in 29] is discovered a so-called matrix phase, which is not present in the continuous 
model, as a result the infinite matrix limit fails to converge to the classical case. The adding 
of an extra term to the matrix action |3(| solves the problem making the new anomalous phase 
transition to disappear. The extension of the parameters space, together with the classification 
of the different phase regions, would allow us to compare our model with the results of the 
simulation performed for the fuzzy sphere, have a look at the occurrence of a correspondent 
matrix phase and eventually to try to regularize the infinite matrix limit. 
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